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(^-HARMONIC  FUNCTIONS  AND  INVERSE  CONDUCTIVITY  PROBLEMS 

ON  NETWORKS. 


CARLOS  A.  BERENSTEIN  AND  SOON-YEONG  CHUNG 


Abstract.  In  this  paper,  we  discuss  the  inverse  problem  of  identifying  the  connectivity  and 
the  conductivity  of  the  links  between  adjacent  pair  of  nodes  in  a  network,  in  terms  of  an  input- 
output  map.  To  do  this  we  introduce  an  elliptic  operator  and  an  cu-haxmonic  function  on  the 
graph,  with  its  physical  interpretation  been  the  diffusion  equation  on  the  graph,  which  models 
an  electric  network.  After  deriving  the  basic  properties  of  u;-harmonic  functions;  we  prove  the 
solvability  of  (direct)  problems  such  as  the  Dirichlet  and  Neumann  boundary  value  problems. 
Our  main  result  is  the  global  uniqueness  of  the  inverse  conductivity  problem  for  a  network  under 
a  suitable  monotonicity  condition. 


Introduction 

A  network  represents  a  way  of  interconnecting  any  pair  of  users  or  nodes  by  means  of  some 
meaningful  links.  Thus,  it  is  quite  natural  that  its  structure  can  be  represented,  at  least  in  a 
simplified  form,  by  a  connected  graph  whose  vertices  represent  nodes  and  whose  edges  represent 
their  finks.  When  we  have  some  problem  on  a  part  of  the  network  or  when  we  are  in  need  of 
finding  such  problem,  it  is  almost  impossible  to  investigate  the  whole  network,  since  the  network 
may  be  too  vast  and  its  structure  or  connectivity  too  complicated. 

From  the  graph  theoretical  point  of  view,  problems  involving  graph  identification  have  been 
among  the  most  important  and  famous  open  problems  in  graph  theory  (  [BH]  ).  Most  of  the 
work  on  this  subject  has  concentrated  on  spectral  graph  theory,  on  the  realization  of  graphs  with 
given  distances,  and  on  the  reconstruction  of  graphs  from  vertex  deleted  subgraphs  (  see  [B2], 
[Cl].  [C2],  [C3],  [CO].  [CL],  [CvDGT],  [CvDS]  and  [HY]  ).  Thus  far.  spectral  theory  has  been 
one  of  the  most  significant  tools  used  studying  graphs,  and  it  has  led  to  noteworthy  progress  in 
the  study  along  these  questions.  But,  as  it  is  well  known,  graphs  are  not  in  general  completely 
characterized  by  their  spectra  (see  [CvDGT],  p.  66). 

In  this  paper  another  method  to  study  the  graph  identification  problem  will  be  introduced,  a 
discrete  version  of  the  inverse  conductivity  problem. 

The  inverse  conductivity  problem  original  aim  was  to  identify  the  conductivity  coefficient  in 
continuous  media  from  boundary  measurements,  such  as  Dirichlet  data,  Neumann  data,  or  their 
appropriate  combinations. 
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The  discrete  or  finite  nature  of  graphs  makes  working  on  graphs  basically  easier  than  inves¬ 
tigating  these  problems  in  the  continuous  case.  On  the  other  hand,  their  discrete  nature  of  also 
gives  rise  to  several  disadvantages.  For  example,  solutions  of  the  Laplace  equation  (introduced  in 
Section  2)  have  neither  the  local  uniqueness  property  nor  is  their  uniqueness  guaranteed  by  the 
Cauchy  data,  contrary  to  the  continuous  case,  where  they  are  the  most  important  mathematical 
tools  used  to  study  the  inverse  conductivity  problem  and  related  problems. 

The  purpose  of  this  paper  is  to  give  a  discrete  analogue  of  the  inverse  conductivity  problem 
studied  as  studied  in  a  number  of  publications,  such  as  [A].  [Ca],  [I],  [IP],  [KS]and  [SU].  To  do 
this  we  introduce  an  elliptic  operator  on  the  graph,  the  uj  -  Laplacian  and  interpret  it  as 
a  diffusion  equation  on  the  graph  modeled  by  the  electric  network.  Since  little  has  been  so  far 
studied  about  partial  differential  equations  on  graphs,  we  will  establish  several  useful  properties 
of  Aw,  which  are  essential  to  solve  the  inverse  problem. 

The  inverse  problem  we  study  is  to  identify  the  connectivity  of  the  nodes  and  the  conductivity 
of  the  edges  between  each  adjacent  pair  of  nodes.  We  begin  by  proving  the  following  global 
uniqueness  result  for  the  inverse  conductivity  problem  in  a  network  satisfying  the  monotonicity 
condition: 


Theorem.  Let  and  002  be  weights  with  ooi  <  002  on  S  x  S  and  /1,  f‘2  :  S  — >  M  be  functions 
satisfying  that  for  j  —  1,  2, 

A^/jO)  =  0.  x  G  S, 

<  grW)=UA  *  ^dS; 

.  Is  'fj<Li=K 

for  a  given  function  f)  :  dS  — >  M  with  Jds  =  0  and  for  a  suitably  chosen  number  K  >  0. 

If  we  assume  that 

o 

(i)  uj\{z.y)  =  uj2{z.y)  on  dS  x  dS. 

(ii)  /i|as  =  /2|as  ? 

then  we  have 

fi  =  h  on  S 


and 


uji  =  0J2  on*?  x  S. 


The  second  conclusion  uj\  =  u>2  above  is  exactly  what  we  want  to  have.  In  fact,  it  shows  not 
only  whether  or  not,  each  pair  of  two  nodes  is  connected  by  a  link,  but  also  how  nice  the  link  is. 

Both  the  condition  u)\  <  tV2  above  (so  called,  the  monotonicity  condition  )  and  the  condition 
fs  fj  d^  —  K  (so  called  the  normalization  condition)  will  be  shown  to  be  essential  by  giving 
counterexamples.  In  fact,  even  in  the  continuous  case,  some  form  of  monotonicity  has  also  been 
considered  (  see  [I],  [Ca]  and  [A]  ). 

We  organized  this  paper  as  follows:  First,  we  discuss  calculus  on  graphs  in  Section  1  and  in 
Section  2  we  introduce  cj-harmonic  functions  on  graphs  and  some  good  properties  of  them,  which 
are  useful  later  and  for  further  study.  In  fact,  those  properties  are  interesting  by  themselves  in 
authors*  opinion. 


w-HARMONIC  FUNCTIONS  AND  INVERSE  CONDUCTIVITY  PROBLEMS  ON  NETWORKS. 


3 


In  Section  3.  we  discuss  the  direct  problems  such  as  the  Dirichlet  BVP  and  Neumann  BVP. 
and  give  a  physical  interpretation  of  Aw.  Besides,  additional  useful  properties  of  ^-harmonic 
functions  will  be  introduced. 

In  the  final  Section  4.  we  prove  the  global  uniqueness  result  of  inverse  problem  under  the 
monotonicity  condition.  Ahead  of  its  proof,  we  derive  an  discrete  version  of  the  Dirichlet  prin¬ 
ciple.  which  is  an  essential  tool  for  the  proof  of  the  main  theorem. 

After  the  authors  completed  this  paper.  Professor  Gunter  Uhlmann  informed  the  authors  that 
Morrow  with  his  group  published  a  series  of  papers  (  see  [MCI].  [MC2].[MMC].  [MI]  and  [MIC] 
)  on  the  inverse  problem  of  the  networks.  But  their  results  were  concentrated  on  the  networks 
of  special  types  such  as  circular  networks  or  integer  lattices.  Moreover,  their  approaches  would 
not  work  for  the  networks  of  general  type. 

1.  Calculus  on  Weighted  Graphs 

We  shall  begin  with  some  definitions  of  graph  theoretic  notions  frequently  used  throughout 
this  paper. 

By  a  graph  G  =  G{V.E)  we  mean  a  finite  set  V  of  vertices  with  a  set  E  of  two-element 
subsets  of  V  (whose  elements  are  called  edges).  The  set  of  vertices  and  edges  of  a  graph  G  are 
some  times  denoted  by  V{G)  and  E{G).  or  simply  V  and  E.  respectively.  But  conventionally, 
we  denote  either  x  G  V  or  x  G  G  the  fact  that  x  is  a  vertex  in  G. 

A  graph  G  is  said  to  be  simple  if  it  has  neither  multiple  edge  nor  loop  and  G  is  said  to  be 
connected  if  for  every  pair  of  vertices  x  and  y  there  exist  a  sequence  (termed  a  path)  of  vertices 
x  =  to.ti.T2-***  * xn~\ . xn  =  y  such  that  Xj-i  and  xj  are  connected  by  an  edge  (termed 
adjacent)  for  j  =  1. 2.  •  •  •  .  n. 

A  graph  S  =  S{Vf.  Er)  is  said  to  be  a  subgraph  of  G(V.  E)  if  Vf  C  V  and  Er  C  E.  In  this  case, 
we  call  G  a  host  graph  of  S.  If  Ef  consists  of  all  the  edges  from  E  which  connect  the  vertices 
in  Vf .  then  S  is  called  an  induced  subgraph.  It  is  easy  to  see  that  every  induced  subgraph  of  a 
connected  graph  is  also  connected. 

A  weighted  (undirected)  graph  is  a  graph  G{V.  E)  associated  with  a  weight  function  uj  : 
V  x  V  — >  [0.  oc)  satisfying 

(i)  uj(x.x)  =  0.  x  G  V. 

(ii)  uj(x.  y)  =  oj{y*  t).  if  x  ~  y. 

(iii)  t o{x.y)  =0  if  and  only  if  {x.  y}  i  E. 

Here,  x  ~  y  means  that  two  vertices  x  and  y  are  connected  (adjacent)  by  an  edge  in  E.  In 
the  case.  {x.  y}  denotes  the  edge  connecting  the  vertices  x.  y. 

In  particular,  a  weight  function  uj  satisfying 

uj{x.y)  =  1.  if  x  ~  y 

is  called  the  standard  weight  on  G.  The  physical  meaning  of  the  weight  function  will  be  discussed 
later  in  Section  3. 

The  degree  d^x  of  a  vertex  x  in  a  weighted  graph  G{V .  E)  with  a  weight  uj  is  defined  to  be 

(Lx  :=  yYO ,y). 
yev 
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Throughout  this  paper,  all  the  subgraphs  are  assumed  to  be  induced  subgraphs  of  a  host 
graph,  which  is  simple  and  connected,  with  a  weight,  and  a  function  on  a  graph  is  understood 
as  a  function  defined  only  on  the  set  of  vertices. 

The  integration  of  a  function  /  :  G  — >  R  on  a  graph  G  —  G(V>  E)  is  defined  by 


(or  simply  /  /)  :=  Y  f{x)dUJx. 
Jg  xev 


We  shall  now  define  the  directional  derivative  of  a  function  /  :  G  — ►  M.  For  each  x  and  y  £  V 
we  define 

~im  . 

The  gradient  of  function  /  is  defined  to  be  a  vector 


Vw/(ar)  :=  (D^yf(x))yeV  , 

which  is  indexed  by  the  vertices  y  £  V.  Then  it  is  easy  to  see  that 


iv^/(x)i2^ 

E  E  I Ay)  -  /0)  I2  u{x,  y ) 

v  yev 

2  E  I  f(v)~  f(x)G(x,y), 

{x,y}£E 


which  is  called  the  energy  of  /  on  G . 

For  a  subgraph  S  of  a  graph  G  =  G{V.  E)  the  (vertex)  boundary  dS  of  S  to  be  set  of  all 
vertices  z  G  V  not  in  S  but  adjacent  to  some  vertex  in  S .  i.e. 


dS  {z  G  V\z  y  for  some  y  G  S} 

o 

and  we  define  the  inner  boundary  dS  is  defined  by 

o 

dS:=  {y  G  S\y  for  some  z  G  dS} . 
Also,  by  S  we  denote  a  graph  whose  vertices  and  edges  are  in  S  U  dS. 


The  (outward)  normal  derivative  ^j^(z)  at  z  G  dS  is  defined  to  be 


#C>  :=  Ei«2)  - /w  ■ 

U0J  1 1  e 


y(z,y) 

dLz 


where  d'^z  =  'Eyes  u’(z-  ]))■ 

The  to-Laplacian  of  a  function  /  :  G  — >  M  on  a  graph  G  is  defined  to  be 


A uf(x)  :=  -  Yi  ( Dlf )  (x).  V. 

yev 
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In  other  words. 

(1-1)  A^/0)  :=  J^[/(y)  -  f(x)]  ■  ,  x  G  V. 

yev 

For  notations,  notions  and  conventions  we  refer  to  [Cl]  and  [CvDS]. 

Remark  1.1.  (i)  The  discrete  Laplacian  on  graphs  can  be  found  in  several  places,  such  as 

[Cl].  [CvDS].  [Bl].  But  the  a;-Laplacian  defined  above  is  not  exactly  the  same  as  the  one 
considered  in  those  references.  In  fact,  the  definition  used  here  will  give  us  an  advantage 
of  a  more  consistent  treatment  in  Section  2. 

(ii)  The  first  derivatives  and  gradient  in  a  discrete  sense  have  not  been  introduced  so  far 
precisely  in  the  literature,  as  far  as  the  authors  know.  But  the  first  derivative  D^.y 
defined  above  may  still  be  unsatisfactory  in  a  sense  that  Leibniz*  rule  does  not  hold.  In 
spite  of  this  defect,  it  will  be  seen  later  that  it  has  the  appropriate  physical  meaning  and 
works  very  well  with  respect  to  calculus  on  graphs. 

In  what  follows,  a  function  /  defined  on  S  may  be  understood  as  a  function  on  its  host  graph 
G  such  that  f  =  0  on  G\SAi  necessary. 

Theorem  1.2.  Let  S  be  a  subgraph  of  a  host  graph  G.  Then  for  any  pair  of  functions 
f  :  S  ^  M  and  h  :  S  — >  M,  we  have 

(1.2)  2  [_  hi-A^f)  =  [  V^h  •  Vw/. 

Js  Js 

Proof.  A  direct  use  of  the  definitions  mentioned  above  gives 

2  f_ M-AV)  =  T^)[-A VO)K^ 

&  x(zS 

=  ~^Yh^{  Y  f(x)Mx-y)} 

xeS  yeV(G ) 

=  2YYh(-x^^~  /(y)M^y) 

xES  yES 

=  Y  Y  h(x)[f(x)  -  f(y)Mx->  y)  +  YY  h(v)[f(v )  -  /0)NY  y) 

xES  yES  xES  yES 

=  YY  { [Ay)  -  /(*)]  }  •  { [Ky)  -  Kx)]  a/vAV } 

xES  yES 

=  Y  {W/0)  ■  Vojh^jdajX 

xES 

=  jfvw/-V<A  ■ 

Js 

The  above  theorem  yields  many  useful  formulas  such  as  the  Green  theorem. 

Corollary  1.3.  Under  the  same  hypotheses  as  above  we  have  the  following  identities: 
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(i) 


2  £/(- Awf)=  [  |W/|2. 

Js  Js 


(ii) 


[hA^f  =  f  fA^h. 

Js  Js 


(iii)  (Green’s  formula) 

[  (fAU)h  -  hAJ-)  =  f  (/|^  -  hp-). 

Js  JdS  dojU  djn 

Proof,  (i)  is  trivial  and  (ii)  can  be  easily  obtained  by  the  symmetry  in  (1.2).  We  prove  (iii).  In 
view  of  (ii)  we  have 

0  =  [  [fAJi,  -  hAuf] 

Js 

=  [  [fAcjh  -  hAuf]  +  [  [fA^h  -  hAuf] 

Js  JdS 

Then,  since  S  is  the  induced  subgraph,  it  follows  that  uj{z.  y)  =  0  for  all  z  and  y  £  dS  and 

({fA^h-hA^n  -  [  [hA^f-fA^h] 

Js  JdS 

=  ^2  -  f(z)A0Jh(z)]d0Jz 

zeds 

51  {h(d  [. f(y )  -  /(*)]  w(z.  y)  -  f(z)]  \h(y)  -  h(z)\  w(z.  y) } 


zeds  yes 


51  [h{z){ 


zeds 


df  \  ,  tf  \  dh  f  \ 

a^}+/(2Wz) 


d^z 


=  L 

JdS  Quid'  ojn 


In  the  continuous  case,  the  followings  are  well-known  formula  : 

A{  fg)  =  fAg  +  2 V/  •  Vg  +  gAf 

[  Vf-Vg+  [  fAg  =  f  f^ 

Jn  Jn  Jdn  ^ 

Here,  we  introduce  a  discrete  analogue  of  these  formula. 


Theorem  1.4.  Under  the  same  hypothesis  as  in  Theorem  1.2,  we  have  following  identities  hold 
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(i) 


At —  f Awh  +  Vw/  •  Vu/i  +  hA^f 


(ii) 

f  Vu,f-Vu,h  +  [  IfAJi  +  hAJ]  =  [ 

Js  Js  JdS  <h»- 

Proof,  (i)  can  be  obtained  by  an  elementary  manipulation.  Using  now  (i)  and  Theorem  1.2.  (iii) 
with  h  =  1  we  obtain  (ii).  g 


2.  uj  -  Harmonic  Functions 

In  this  section  we  will  discuss  the  functional  properties  of  functions  which  satisfy  the  equation 

(2-1)  A uf(x)  :=  Y  [ f(y )  -  /(*)]  -  =  °- 

yes 

For  a  subgraph  S  with  boundary  dS  ^  <j>  of  a  host  graph  G  with  a  weight  uj  we  say  that  a 
function  /  :  S'  — ►  R  is  uj-harmonic  on  S'  if  it  satisfies  (2.1)  for  all  x  £  S.  i.e. 

/V)  =  f—  Y  y)-  x^s. 

yes 

This  implies  that  the  value  of  /  at  x  is  given  by  a  weighted  average  of  the  values  of  /  at  its 
neighboring  vertices.  From  this  point  of  view,  we  can  clearly  expect  the  following  result  to  be 
true: 

Theorem  2.1  (Minimum  and  Maximum  Principle).  Let  S  be  a  subgraph  of  a  host  graph  G  with 
a  weight  uj  and  f  :  S  ^  M  be  a  function. 

(i)  If  A  ojf(x)  >  0  .x  G  S  and  f  has  a  maximum  at  a  vertex  in  S}  then  f  is  constant. 

(ii)  IfAufix)  <  0.  x  £  S  and  f  has  a  minimum  at  a  vertex  in  S,  then  f  is  constant. 

(iii)  If  A0Jf(x)  =  0.  x  £  S  and  f  has  either  a  minimum  or  maximum  in  S}  then  f  is  constant. 

(iv)  If  A(jJf(x)  =  0.  x  £  S  and  f  is  constant  on  the  boundary  dS ,  then  f  is  constant. 

Proof,  (ii)  can  be  done  in  a  similar  way  as  in  (i).  (iii)  and  (iv)  are  easily  obtained  from  (i)  and 
(ii). 

We  prove  (i).  Assume  that  /  has  a  maximum  at  a  vertex  x$  £  S.  Then 


(2.2) 

/Vo)  >  f(y),  y^S 

and 

(2.3) 

yes 
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Suppose  that  there  exists  yo  G  S  such  that  ~  y$  and  f(x o)  ^  /(yo),  i.e.  f(x o)  >  f(yo)  in 
view  of  (2.2).  Then  it  follows  from  (2.3)  that 


f(x  o)  < 


E 

yes 
y^y  o 


f{y)to{x^y)  f(yo)Lo(x0,yo) 

dojXo  d^xo 


< 


E 

yes 

y^yo 

=  /O  o)? 


f(xo)w(x0,y)  ,  f(x0)w(xo,yo) 


d^x  o 


+ 


d^XQ 


which  implies  that  f(x o)  =  f(y)  for  all  y  £  S'  such  that  y  ~  xq.  Now  for  any  x  £  S',  there  exists 
a  path 

xo  ~  xi  ^  x 2  ^  *  *  *  ^  xn—i  ^  xn  —  x. 

since  S  is  connected.  By  the  applying  the  same  argument  as  above  inductively  we  see  that 
/  Oo)  =  f(x).  ■ 


The  following  is  an  easy  consequence  of  the  above  theorem 

Corollary  2.2.  Under  the  same  hypotheses  as  in  Theorem  2.1,  the  following  statements  are 
true: 

(i)  If  A ^f  >  0  on  S  and  f\$s  <  0  (<  0),  then  f  <  0  (<  0)  on  S. 

(ii)  If  A ^f  <  0  on  S  and  f\os  >  0  (>  0),  then  f  >  0  (>  0)  on  S. 

Corollary  2.3.  (1)  If  two  functions  f  and  g  on  S  satisfy 

A ^f  =  0  and  A^g  >  0 

on  S ,  then  g\os  <  f\dS  implies  y  <  /  on  S. 

(2)  If  a  function  f  :  S  —>  R  satisfies 

A ojf(x)  =0.  x  £  S 

and  |/|  has  a  maximum  in  S ,  then  f  is  constant. 


In  the  continuous  case,  it  is  well  known  that  a  local  maximum  principle  holds  for  a  harmonic 
function  in  an  open  subset  Q  C  MT.  But  it  is  not  hard  to  see  that  the  local  maximum  principle 
is  no  longer  true  in  general  in  our  case.  Moreover,  the  local  uniqueness  principle  does  not  hold 
in  general.  As  a  matter  of  fact,  it  is  rather  natural  to  expect  that  such  discrepancies  are  caused 
by  the  discrete  nature  of  graphs. 

A  nonempty  subset  T  of  vertices  of  a  subgraph  S  is  said  to  be  a  surface  in  S  if  T  =  dT  for  a 

O  O 

subgraph  T  whose  vertices  belong  to  S.  In  this  case,  we  denote  by  T  the  inner  boundary  dT. 

o 

For  each  vertex  z  £  T  and  x  £T  we  define 

diz  :=  ^  to{y.  z)  (inward  degree) 
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and 


d!fx  :=  ^  uj(x.  z)  (outward  degree). 
zer 


In  addition,  for  a  function  f  on  S  we  write 


Y  f(z)diz 

zer 


(inward  integral) 


and 


f(x)c&x 


Y  f(x)d!^x  (outward  integral) 

G 

^er 


We  use  these  notions  to  obtain  the  following  interesting  properties  of  ^-harmonic  functions. 


Theorem  2.4.  Let  S  be  a  subgraph  of  a  host  graph  with  weight  to  and  let  f  :  S  — >  ®L  Then  f 
is  uj-harmonic  on  S,  i.e.}  for  all  x  £  S, 


(2.4)  A^/Or)  =  0, 

if  and  only  if  for  every  surface  T  in  S 

(2-5)  ^  f(z)dp  =  J  f(y)d'^y 

_  _  O 

Proof.  Let  x  G  S  and  Tx  =  [y  G  S\x  ^  yj  .  Then  Tx  is  a  surface  in  S  and  Tx=  {t}.  Since 

o 

d^x  =  drfx  on  Tx  and  d^z  =  uj(x.  z).  (2.5)  implies 

f(x)d^x  =  Y  f(z)v(x,z). 


which  implies  (2.4)  immediately 

Assume  now  that  (2.4)  holds  and  let  T  be  a  surface  in  S'  such  that  T  =  dT  for  a  subgraph  T  C  S. 
We  use  Green's  formula  (Corollary  1.3.  (ii))  to  obtain 


(2.6) 


0  = 


/  A/ 

Jt 


[dp 
r  djn 
dp 
d^n' 


\  ^  df  f 

= 


zer 


=  YY  [#*)  -  f(y)]v(z>y)- 


zer  ° 
ye  r 


Then  it  follows  that 


Y  Y  / OM^  y)  =  Y  Y  y) 
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or.  equivalently 


which  yields  (2.5). 


^2f(z)  ='52f(y) 


zer 


yer 


ye  r 


zer 


In  view  of  (2.6)  we  obtain  the  edge  version  of  Theorem  2.4.  the  so  called  dual  theorem,  as 
follows  : 


Corollary  2.5.  Under  the  same  conditions  as  in  Theorem  2.4  ,  the  formula  (2.5)  is  equivalent 
to 

Y  [f(z)~  f(y)\u(z,y)  =  0 

{x,y}eE(  r,f) 

o  o 

where  £?(T.  r)  denotes  the  set  of  all  edges  joining  a  vertex  in  T  and  a  vertex  in  T- 

For  two  vertices  x  and  y  in  a  connected  graph,  the  distance  d{x.  y)  between  x  and  y  is  the 
number  of  edges  in  a  shortest  path  joining  x  and  y. 

For  a  vertex  x$  in  a  subgraph  S'  we  write 

UOro)  :=  { V  G  S\  d(x0)  -j}  ,j  =  0, 1,2,  ••• 

which  is  called  a  neighborhood  of  xq  with  radius  j. 

Then  the  following  is  a  variant  of  Theorem  2.4: 

Corollary  2.6.  Let  S  and  f  be  the  same  as  in  Theorem  2.4 ■  Then  f  is  uj-harmonic  on  S  if 
and  only  if  for  every  x$  G  S 

(2.7)  [  f(x)d£x  =  f  f(x)d!wx 

JTj{x  0)  JTj+xixn) 

for  each  j  with  Tj(xo)  C  S. 

Proof.  Letting  j  =  0  in  (2.7)  we  have  the  sufficiency.  To  prove  the  necessity,  consider  an  induced 
subgraph  T  whose  vertices  are  exactly  those  of  Ui=o  ^(^o)-  Then  it  is  easy  to  see  that 

dT  =  rj+i(xo)  and  dTc  Tj(x o). 

o 

But  a  vertex  x  in  Tj(x o).  which  does  not  belong  to  dT.  does  not  make  any  contribution  to  the 
outer  integral  /r.^0^  f(x)d^x.  since  drfx  =  0.  Hence,  condition  (2.5)  in  Theorem  2.4  shows  the 
condition  is  necessary,  g 


The  following  is  the  dual  version  of  the  above  corollary: 

Corollary  2.7.  Under  the  same  conditions  as  in  Corollary  2.6  the  formula  (2.7)  is  equivalent 
to  _ 

Y  l/U)  -  f(y)]  y)  =  o 

{x,y}eE  (r  j  (a?o).r  j+1  (a?o)) 

where  E(Tj(xo).  Tj+i(xo))  denotes  the  set  of  all  edges  joining  a  vertex  in  Tj(xo)  and  a  vertex 
in  ri+i(x0) 
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3.  The  Dirichlet  and  Neumann  Boundary  Value  Problems.  Direct  Problems 


We  start  this  section  with  a  physical  interpretation  of  the  a;-Laplace  and  a;-Poisson  equations. 
Consider  a  host  graph  G  with  a  weight  uj  and  an  (induced)  subgraph  S.  For  a  surface  r  in  S 
with  r  =  dT  for  some  T  C  S  and  z  G  T.  the  flux  of  energy  passing  through  z  to  its  adjacent 
nodes  in  T  is  given  by 

(3-D  -E[ 

yr^z 

where  d!z  =  Yy^yeT^^'V)  and  /  is  a  potential  function  in  a  diffusion  field  on  a  network. 
(For  example,  an  electrostatic  field,  a  thermal  field,  or  an  elastic  membrane.)  Here,  the  weight 
uj(z.y)  plays  the  role  of  the  conductivity  of  the  diffusion  along  the  edge  {z.y}.  In  fact.  (3.1)  is 
exactly  —  J^Oz)  on  r  by  definition  (see  Section  1)  and  thus,  by  Green's  formula  we  have 

which  is  the  flow  across  T. 

On  the  other  hand,  assume  that  T  gains  (or  loses)  an  amount  of  energy  fT  g  where  g  is  the 
energy  density.  Then  we  have 

/  (-A,/)  =  [  9 

JT  JT 

Therefore,  since  T  is  arbitrary,  by  taking  T  to  be  any  single  vertex  x  G  S  we  obtain  the  vertex 
equation 

(3.2)  -A w/(x)  =  g{x),  x  G  S. 


Thus,  it  is  reasonable  to  say  that  the  conductivity  equation  on  a  graph  can  be  represented  as  in 

(3.2).  where  uj{x.y)  corresponds  to  the  edge  conductivity  on  the  edge  x.y. 

Following  the  work  of  Fan  Chung  and  her  collaborators  [Cl].  [C2]  and  [CO],  we  wifi  discuss 
first  the  equation  (3.2)  on  a  graph  G  =  G(V*  E)  with  a  weight  uj  and  no  boundary.  We  consider 
the  matrix 

[  0.  otherwise  . 

Considering  the  function  /  as  a  |V|-dimensional  vector,  the  equation  (3.2)  can  be  understood 
as  a  matrix  linear  equation.  Let  D  denote  the  diagonal  matrix  with  (x.  x)-th  entry  having  the 
value  d^x  for  each  x  and  Cc  =  D^AuD-U*.  Then  (  — Cuj)  is  a  nonnegative  definite  symmetric 
matrix,  so  that  it  has  the  eigenvalues 


Aq  <  Ai  <  A2  <  ■  ■  ■  <  Aw-i 


and  the  corresponding  eigenfunctions 

(3.3)  $0;  ^1;  '  '  '  <&AT— 1* 
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which  form  an  orthonormal  basis  for  R.N  in  the  sense  that  for  each  pair  of  distinct  i  and  j 

Y  *i(ar)  •  $j(x)  =  0, 


while,  for  all  j. 


E  =  i- 


Here.  N  denotes  \V\.  the  number  of  vertices  in  G.  Then  it  is  easy  (see  [Cl])  to  show  that 
A0  =  0,  Ai  >  0  and  3>o0*0  =  ,  x  G  V,  and  vol(G)  :=  ^2xeV 

V  vol(G) 

In  what  follows,  we  occasionally  use  the  notation  (  .  }x*  defined  by  ( f.g)x  —  Yhxex  f(x)9(x) 
for  simplicity.  Now  we  have  the  following  solvability  result  for  the  Poisson  equation: 

Theorem  3.1.  Let  G  =  G{V .  E)  be  a  graph  with  a  weight  uj  and  f  :  G  — ►  R  be  a  function.  Then 
the  equation 

(3.4)  Y/Or)  =  g{x),  x  6  V 

has  a  solution  if  and  only  if  fG  g  =  0.  In  this  case ,  the  solution  is  given  by 

(3.5)  f(x)  =  a0  +  (rw(x,  ■),  g)v,  x  €  V 

where  ao  is  an  arbitrary  constant  and 

N- 1 

(3.6) 


T^{x,y)  =  Y  x.,y£V. 


i= i 


Proof.  Assume  that  fc  g  —  0.  Then 


(■ D1/2g ,  $0)  =  Y  VXxg(x) 


X 


=  f 

volG  Jg 


\J volG' 


Vvol G  jg 
=  0 


9 


where  D  is  the  diagonal  matrix  whose  x-th  diagonal  entry  is  dwx. 

Consider  the  orthogonal  expansion 

N-l 

(■ D1/2f)(x )  =  Y  X  e  V 

3= 0 

where  aj  =  ( D 1^2/.  &j)*j  =  0. 1. 2.  •  •  •  .  TV  —  1.  Then  since  C^D1^2  =  D1^2 and 

-Ajaj  =  (D^f.C^j) 

=  (CvDWf.Qj) 

=  {D'Pg^j), 


=  (-()(G2Gl  j  =  1,2,---  .N-l 


we  have 
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and  ao  is  an  arbitrary  constant.  Hence 

\fdojX 
a/voIG 


vVA/O  =  Q0  +  9(y)$j(y)V<isy  $j(x) 

j= l  ^  ytV 


equivalently. 


which  gives  (3.6). 

The  proof  of  the  converse  is  easy. 


The  matrix  Tw  in  (3.5)  is  called  the  Green  function  of  Aw. 

The  following  corollary  is  a  Liouville  type  theorem  for  ^-harmonic  functions. 

Corollary  3.2.  Under  the  same  conditions  as  in  Theorem  3.1 ,  every  solution  f  of 

Au>f(x)  =  0.  x  G  V 

is  constant. 


The  following  corollary  describes  all  functions  which  are  ^-harmonic  except  possibly  on  a 
given  (singularity)  set  T. 

Corollary  3.3.  Under  the  same  conditions  as  in  Theorem  3.1 }  let  T  C  V.  Then  every  solution 
to 

A „f{x)  -  0,  x  G  V  \  T 

can  be  represented  as 

(3.7)  f{x)  =  a0  +  y]ra,(x.y)o;(y).  x  G  V 

yeT 

where  ao  is  an  arbitrary  constant  and 

a{y)  =  Aufiy).  yeT. 

In  particular,  ii  T  =  {xq}.  x$  G  V.  then  (3.7)  can  be  written  as 

f(x)  =  a0  +  aoT^ix,  x0).  x  G  V 

where  a0  =  A w/(t0). 

Let  us  now  turn  to  boundary  value  problems  and  their  eigenvalues.  For  a  subgraph  S'  of  a 
host  graph  G  with  a  weight  c o.  the  Dirichlet  eigenvalues  of  —  =  —D1/2 A^D-1/2  are  defined 
to  be  the  eigenvalues 

ZY  <  v2  <  *  *  *  <  isn 
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of  the  matrix  —C^.s  where  C^.s  is  a  submatrix  of  with  rows  and  columns  restricted  to  those 
indexed  by  vertices  in  S'  and  n—  |S|.  Let  0u02c*  ;0n  be  the  functions  on  S  such  that  for 
each  j  =  1,  2,  *  *  *  ,  n, 


£u,S<t>j{x)  =  x  G  S  and  4>j\ds  =  0. 

In  fact.  0i,  02-  ■  *  *  *  0n  are  the  eigenfunctions  corresponding  to  z^i  <  z/2  <  *  *  *  <  and  form  an 
orthonormal  basis  for  ML.  Then  it  is  easy  to  verify  that  the  first  eigenvalue  v\  >  0,  (see  for 
instance,  [Cl]). 

One  can  follow  now  the  standard  procedure  to  define  Green  functions  7 ^.5  as  follows  : 


(3.8) 


|S| 

lu,s(x,y)  =  Y  ( 
3= 1 


V3 


Vdu/y 

\fdjc' 


x.y  G  S 


A 


Letting  stand  for  the  diagonal  matrix  whose  x-th  entry  is  for  each  x  G  S  and  setting 
—  Ds^CvsD]/2.  one  can  easily  verify  that 


(3.9)  'Yu>,s\j,s  =  \j,slu>,s  =  I 
and 

151  A/ 

(3.10)  A^sOr.y)  =  Y=,  x,y  €  S. 

j= 1  V  ^ 

where  /  denotes  the  |S|-dimensional  identity  matrix. 

The  Dirichlet  boundary  value  problem  was  solved  by  F.R.K. Chung  in  [C2],  when  the  graph 
has  the  standard  weight.  (For  the  interested  reader,  despite  some  minor  errota,  the  proof  given 
there  is  correct.)  We  prove  now  the  solvability  of  the  Dirichlet  boundary  value  problem  for 
graphs  with  arbitrary  weights  using  a  different  method. 

Theorem  3.4.  Let  S  be  a  subgraph  of  a  host  graph  with  a  weight  uj  and  a  :  dS  ^  M  be  a  given 
function.  Then  the  unique  solution  f  to  the  Dirichlet  boundary  value  problem  (DBVP) 

f  A uf(x)  =  0.  x  €  S. 

(  /las  =  o’ 


can  be  represented  as 

(3.11) 

where 


f{x)  =  -(7 v(x,  ■),  Ba)yeS,  x  G  S, 


Ba(y )  = 


E 

zEdS 


<r(z)u(y-,z) 

du>y 


(3.12) 


yGS 
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Proof.  Let  /  be  a  solution  of  DBVP.  Then 

(3.13)  0  =  ^  7 UJ.s(x,y)AUJf(y) 

yes 


Here,  we  have  used  Green's  formula  from  the  Corollary  1.3.  On  the  other  hand,  one  can  show 
that 

M^1/2<^)U>  =  (-^)(T?1/2A)(A  x  e  s> 

since  4>j  =  0  on  dS  and 

\s\  nr~ 

^(f>j{x)(t>j{y)J-pr=5{x,y),  x.  y  €  S. 

■  1  V  CL(jjX 

J= 1 

where  S  denotes  the  Kronecker  delta.  From  these  identities  we  can  conclude  that 
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Hence,  from  the  equality  (3.13)  and  the  fact  that  —  0  on  dS.  we  have 

\s\  1  ,,  V 

1  Mx) 


dtjn v  D  ' 


m  =  E(-^P  / 

|S|  1  ^  - 


A'  isj  WdujX 


zEdS 


|S| 


i  Mx) 


E(--)  ^ 

j= 1  3  Vduj  zeds 


E  a(z)dz[J2i 


Mz )  A(y)|wGy) 


yes 

E(C)$gE^)^(E 


|S| 


GC?  GGr  ^  -l 

a{z)io{z,y) 


zEdS 


cLy 


-) 


J=1  J  v  ygg 

=  -^iu.s{x-.y)B<r(y) 

tjES 

=  -(l^.s(X‘-)-,Bcr)s 

for  each  x  £  5. 

The  desired  uniqueness  result  now  follows  easily  from  Theorem  2.1. 


Remark  3.5.  (i)  The  identity  (3.11)  can  be  rewritten  as 


<r{z)u(y,z)- 

(Ly 


4>j(y)<t>jix) 


x  G  S. 


In  fact.  Bc r  is  a  function  on  S  depending  only  on  the  value  of  a  on  dS  and  Ba(y)  =  0  for 

O 

y  £  S\  dS.  On  the  other  hand,  two  different  boundary  conditions  <j\  and  <12  may  give 
rise  to  the  same  solution  whenever  Bai  =  Ba2 . 

(ii)  (3.11)  can  be  understood  as  a  matrix  multiplication  by 

(3.14)  /  =  -  Ju>,s  ■  Ba  on  S 
or.  equivalently. 

(3.15)  A u.sf  =  ~Ba  on  S' 

in  view  of  (3.8).  The  relation(3.15)  enables  us  to  identify  uniquely  the  boundary  values 
from  a  oi-harmonic  function  /  with  A =  0  on  S. 


Now  we  characterize  the  ^-harmonic  functions  with  a  set  of  singularities  in  a  subgraph  with 
nonempty  boundary. 

Theorem  3.6.  Let  S  be  a  subgraph  of  a  graph  with  weight  tv  and  T  C  S'.  Then  every  f  :  S  — ►  C 
satisfying 

Au>f(x)  =  0.  x  G  S  \  T 
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can  be  uniquely  represented  as 

(3.16)  f{x)  =  h{x)  +  Y  7 u>,six- y)P(y)‘.  x  e  s- 

yeT 

where  h  is  a  uj-harmonic  function  on  S  satisfying  h\$s  —  f\dS  and  P(v)  —  AVGA  y^T. 
Proof.  The  uniqueness  is  easy,  by  Theorem  2.1.  Now  let  f3(y)  :=  A0Jf(y).y  £  T.  Then  we  have 


A  ujf{x)  = 


0,  x€S\T, 
f3{x).  igT. 


Define,  for  x  £  S. 


hix)  ■■=Yj^.s(x,y)P(y), 

yeT 


and 


h{x)  :=  fix) -hix). 


Then  h\gs  =  f\dS  and  for  each  x  G  S. 


|S| 


A ujhix)  =  A ufix)  -  Aaj  [  Y  Y  (  -  ~ 

yeT  j= l  ZA 


f(h=  ■  (t>jiy)\fdjypiy) 


|S| 


<W 

,T 


yeTj= i 

Ao,/(x)  -  Y  S(x-  V)0(v) 


4>j(y)V(Ly(3(y) 


yeT 


=  o, 

which  completes  the  proof. 


Remark  3.7.  (i)  In  particular,  if  T  —  {to}  .  To  £  S.  then  (3.16)  can  be  written  simply  as 

fix)  =  hix)  +^.six,  x0)Pix0), 

where  Pix0)  =  A^fix 0). 

(ii)  In  fact,  in  view  of  (3.16)  and  Theorem  3.4.  the  solution  to  the  nonhomogeneous  DBVP 

/  A ofix)  =  gix ).  x  <E  S, 

1  /las  =  a- 

can  be  represented  by 

fix)  =  -{^.six,  ■),  Ba)s  +  Y,s  Y  •)■  g)s- 


Now  we  will  discuss  the  Neumann  boundary  value  problem  (NVBP).  First,  we  recall  Green’s 
formula 


f 


f  9f 
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Hence,  if  there  exists  a  solution  to 

(  A w/  =  g  on  S', 

l  £n  =  ^  °n  dS> 

then  by  Green's  formula  it  is  necessary  that  Jsg  =  Jds  fn 

Theorem  3.8.  Let  S  be  a  subgraph  of  a  host  graph  G  with  a  weight  uj  and  let  f  :  S  — ►  R, 
g  :  S  — ►  R,  and  ijj  :  dS  — ►  R  be  functions  with  fds  f)  —  fsg .  Then  the  solution  to  the  NBVP 

f  A wf(x)  =  g{x),  x  e  S, 

l  £V)=UA  zeds 

is  given  by 


f  (x)  =ao  +  (R(a:,  ■),  g)s  ~  (D(x,  ■),  Z)dS, 

where  T ^  is  the  Green  function  of  on  the  graph  S  as  a  new  host  graph  of  S  and  ao  is  an 
arbitrary  constant. 

Proof.  We  rewrite  (NBVP)  as 

,,  1 f Eyes  [f(v)  -  /(*)]  “g#  =  g(x),  x  G  S, 

[Eyes  [f(v)  ~  fO]  =  “VIA  z  e  9S. 

To  solve  the  system(3.17).  consider  S'  as  a  new  host  graph  with  the  weight  uj  and  with  no 
boundary.  Then  S  is  still  a  subgraph  of  S.  (In  fact,  we  should  note  here  that  if  we  regard  S 

as  a  subgraph  of  G .  then  its  boundary  dS  may  not  be  empty.)  Then,  for  each  z  £  dS.  the 

inner  degree  d^z  is  equal  to  d^z  in  this  new  graph  S.  since  the  induced  subgraph  has  no  edges 
between  the  vertices  on  dS.  Hence  the  equation  (3.17)  can  be  written  as 


(3.18) 


Eyev&  [f(y)  -  f(x )]  =  sOA  xes, 

Eyey0  [f(v)  -  /O]  =  -UA  *  G  dS 


where  Vq  is  the  set  of  vertices  in  S.  Hence  (3.18)  is  equivalent  to 


(3.19) 

where 


V  [f(v)  -  /(A 

yev0 


Lvix^y) 

cZx 


\&(A 


^(cc)  = 


g{x).  x  €  S. 
x  €  dS. 


Therefore,  (NBVP)  is  equivalent  to 

Aw/(a:)  =  ^(x),  x  €  S. 


x  G  S, 
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Thus,  it  follows  from  Theorem  3.1  that 

/0)  =  ao  +  (DO,  •)■  'JCeVo 

=  a0  +  DUy^Q/) 
yev0 

=  ao  +  N  DO.  y)g(y )  -  51  DO:  ^VO) 

yES  zEdS 

=  a0  +  (D0:-):#)S  -  (DO:  •)  V)9S: 
where  ao  is  an  arbitrary  constant.  This  completes  the  proof.  g 

Remark  3.9.  The  solution  to  (NBVP)  is  uniquely  determined  by  the  Neumann  data  ^  on  dS  up 
to  an  additive  constant.  Thus,  we  get  a  unique  solution  if  we  prescribe  the  value  of  /  at  some 
vertex  in  S  or.  for  example,  if  we  seek  the  solution  /  with  Js  f  =  ( a  given  constant). 


4.  Inverse  Problems 


In  the  previous  section,  we  have  seen  that  for  a  function  function  ?/•  :  9S  — >  R  with  Jds  ^  =  0 
the  Neumann  boundary  value  problem 


(NBVP) 


A ajf(x)  =0.  X  G  5, 


has  a  unique  solution  up  to  an  additive  constant.  Therefore,  the  Dirichlet  data  f\os>  z  £  <95  is 
well-defined  up  to  an  additive  constant. 

In  this  section,  we  will  discuss  the  inverse  conductivity  problem  on  the  network  (graph)  S 
with  nonempty  boundary,  which  consists  in  recovering  the  conductivity  (connectivity  or  weight) 
uj  of  the  graph  by  using,  the  so  called  input-output  map.  for  example  by  using  the  Dirichlet  data 
induced  by  the  Neumann  data  (Neumann-to-Dirichlet  map),  with  one  boundary  measurement. 

In  order  to  deal  with  this  inverse  problem  .  we  need  at  least  to  know  or  be  given  the  boundary 
data  such  as  f(z ).  J^(z)  for  z  £  dS  and  uj  near  the  boundary.  So  it  is  natural  to  assume  that 

f\ OS',  gv a\dS  and  ^\ds  ^ s  are  known  (given  or  measured).  But  even  though  we  are  given  all 

these  data  on  the  boundary,  we  are  not  be  guaranteed,  in  general,  to  be  able  to  identify  the 
conductivity  uj  uniquely.  To  illustrate  this  we  consider  a  graph  S  whose  vertices  are  {1.2.3} 
and  dS  =  {0. 4}  as  follows: 


0 

* 


£ 
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with  the  weight 

w(0. 1)  -  1.  w(0.  fc)  =  0  (fc  =  2.  3.  4). 

and 

w(3.4)  =  l.  cj(fc.4)=0  (fc  =  0.  1.  2). 

Let  /  :  5  — >  M  be  functions  satisfying  A0Jf(k)  =  0.  fc  =  1.  2. 3.  Assume  that 
/( 0)  =  0.  /(I)  =  1.  /(3)  =  3,  /( 4)  =  4,  /(2)  =  (unknown). 

Thus,  since  dS=  {1,3},  the  boundary  data  /las,  fjjn 1 0S  and  w|  c  are  known. 
In  fact. 

|{(0)  =  /(0)-/(l)  =  -l, 

|{(4)  =  /(4)-/(3)  =  l. 

od  ' '' 

The  problem  is  to  determine 

cj(1.2)=x.  cj(2.3 )=y.  cj(1.3)  =  z.  and  /( 2). 


From  A0Jf(k)  =0.  k  =  1.2. 3.  we  have 


/(i)  = 

/(2)  = 


=  1. 


/(0)  +  x/(  2)  +  3z 
1  +  r  +  z 

xf(  1)  +  yf{  3) 


x  +  y 

U(i)  +  y/(2)  +  /(4) 


z  +y  + 1 


=  3. 


(4.1) 


/(3)  = 

This  system  is  equivalent  to 

f  x(y  -  1)  +  y{x  -  1)  +  2z(x  +  y)=  0, 


,/(2)  =  31  • 


This  system  has  infinitely  many  solutions.  For  instance,  assume  z  =  0*  that  is.  the  two  vertices 
1  and  3  are  not  adjacent.  Then  (4.1)  is  reduced  to 


(4.2) 


a  + 1  =  2. 

nr  1  ?  /  * 

m  =  gf- 


It  is  easy  to  see  that  there  are  infinitely  many  pairs  (x.  y)  of  nonnegative  numbers  satisfying 
the  first  equation  in  (4.2).  so  that  /( 2)  is  undetermined  as  a  result. 


In  view  of  the  above  example,  in  order  to  determine  the  weight  tv  uniquely  we  need  some  more 
information  than  just  f\ds-,  'n\dS  and  To  m°tivate  the  main  theorem  we  impose  in 

this  example  the  additional  constraints  that 
(4.3)  x  >  1.  y  >  1  and  z  >  0 
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in  (4.1).  Then  the  equation  (4.1)  yields  a  unique  triple  of  solution  x  =  1.  y  —  1.  z  —  0  and 

/( 2)  =  2. 

As  a  matter  of  fact,  even  the  inverse  conductivity  problem  of  a  diffusion  equation  in  a  bounded 
open  subset 


(4.4) 


P[a:u] 


div[a(x)Vu(x)]  =0.  x  G  ft. 
u\dn  =  o- 


has  been  studied  under  some  additional  constraints  besides  Dirichlet  and  Neumann  data  (see 
[A],  [BF],  [Ca].  [I].  [IP],  and  [SU]  ).  In  particular,  in  [A]  and  [I]  it  is  shown  that  there  is  a  global 
uniqueness  result  under  the  condition  that 

(i)  a\  =  <22  near  dQ.  and  a\  <  <22  in  Q. 

(«)  Sf  =  «"  «i. 

(iii)  /,,  »l  („  u2  '-I- 

where  P[aj'.Uj\  =0.  j  —  1.2  in  (4.4). 


Now  we  are  in  a  position  to  state  the  first  main  theorem  of  this  paper. 


Theorem  4.1.  Let  uj\  and  x>2  be  weights  with  uj\  <UJ2  on  SxS  and  f\.  f2  :  S  — >  R  be  functions 
satisfying  that 

(  Au>.fj(x)=0,  x  €  S, 

{  ifn(Z)=^Z)-<  Z^dS 

for  a  given  function  if  :  dS  — >  R  with  f/jS  ip  —  0  and  j  =  1.  2. 

If  we  assume  that 

o 

(i)  toi(z.y)  =u>2(z.y)  on  dS  x  dS. 

(ii)  fl\dS  =  f2\dS  , 
then  we  have 

(i)  fi  =  fi  on  S, 

(ii)  x>i{x.y)  =  x>2(x,y)  whenever  fi(x)  /  fi{y),  or  f2{x)  ^  f2(y). 


To  prove  this  result  we  adapt  the  method  of  energy  functionals,  extensively  used  for  nonlinear 
partial  differential  equations.  For  a  function  a  :  dS  — >  R  we  define  a  functional  by 

(4.5)  Iu[h]:=  f_  [hvuh\ 2-hg\ 

J  s 

for  every  function  h  in  the  set 

(4.6)  A  :=  {/i :  S  — ►  R|  h\os  =  ^ 

which  is  called  the  admissible  set.  In  the  continuous  case,  the  well  known  Dirichlet *s  principle 
states  that  the  energy  minimizer  in  the  admissible  set  is  a  solution  of  the  Dirichlet  boundary 
value  problem.  We  derive  here  the  discrete  version  of  Dirichlet *s  principle  as  follows  : 
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Theorem  4.2  (Dirichlet's  principle).  Assume  that  f  :  S  ^  M  is  a  solution  to 


(4.7) 

Then 


-A  =  g  onS, 

f\ds  =  &■ 


(4.8)  IJ[f]=mmIJ[h}. 

he  A 

Conversely ,  if  f  £  A  satisfies  (4-8),  then  f  is  the  solution  Of  (4-V?  am %  the  only  one. 


Proof.  Let  h  be  a  function  in  A.  Then  making  use  of  (1.2)  in  Theorem  1.2  we  have 

0  =  [(-A^f-gXf-h) 

Js 

=  I  [(-A w/)(/  -  h)  -  g(f  -  h)] 

Js 

=  [_[\vuf-VM-h)-g(f-h)] 

J  s 

=  \  -WA-  j_g{f  -h). 

Hence 


L[\\Vu>f\2  -  9f\  =  L&vf-Vvh-gh] 

J  S  J  s 

<  \  Yi  51 1  ihy)  -  /O]  •  [Ky]  -  h(x)]  I  •  w(ar,  y)~  J  gh 

xeS yeS  ^ 

xES yes  s 

=  iX|v“/|2+iX|v“fc|2-//'*' 

where  we  used  the  triangular  inequality 

.  .  a2  +  h2 

\ab\  <  - *  a?  b  ^ 


Thus.  it  follows  that 


which  implies 
Since  /  £  A.  we  have 


j_[\\V.f\2-gf]<  J_[\\v^h\2-gh\, 

Uf]<Uh]-,  h  €  A. 
min  IJJi]  =  IJ[f}. 
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Xt(x)  = 


Now  we  prove  the  converse.  Let  T  be  a  subset  of  vertices  in  S  and 

1,  xeT 

0.  otherwise. 

Then  /  +  txt  £  A  for  each  real  number  r.  since  xt  =  0  on  dS.  Define 

i(r)  :=  IJ[f  4-  txt],  t  6  M. 

Then 


V)  =  [  [tIVcU  +  tWxtI2  -  (f  +  i~XT)g\ 

J  S 

[  I  Vo,/|2  +  2 rVw/  ■  VvXt  +  t2\VvXt\2  ~  [  (/  +  txt)q 

Js  Js 


1 

4 


Note  that  the  scalar  function  i(r)  has  a  minimum  at  r  =  0  and  thus  ^(0)  =  0.  That  is, 


0  = 


[  W/  -Vo,  xt-  f  xt- g 

Js  Js 


jX 


=  [ Xr(-\jf  -  ^)] 

Js 

=  A./0)  -  g{x)]  gU 

x£T 

In  particular,  taking  T  =  { x }  .  x  G  S.  we  obtain 

-A wf(x)  -  g{x)  =  0, 

which  is  the  required  result.  The  uniqueness  follows  from  Theorem  3.4. 


Now  we  are  ready  to  prove  Theorem  4.1. 

Proof  of  Theorem  4-T 

(i)  Let  a  :  dS  — ►  R  be  the  function  defined  by 

cr(z)  =  fi(z)  =  f2(z),  zeds, 
using  the  hypothesis  (ii).  Define 

V[h]  :=\  [_\V^h\2d^ 

4  Js 

for  every  h  in  the  admissible  set 

A  =  [h  :  S  —>  R\  h\os  =  <r}  • 
Then,  by  virtue  of  Theorem  1.2  we  have 

Ui  [h]  =  1  J  h(—Ai01  h)d^ 

=  \  [  h(-Au}1h)du,1  +  ^  I 

1  Js  1  Jds 
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Moreover,  by  the  coincidence  of  the  Dirichlet  and  Neumann  data  we  can  see  that  the 

o 

boundary  dS  and  the  inner  boundary  dS  are  well-defined  independently  of  the  values  of 
the  weights  uj\.  0J2  and.  moreover,  for  z  G  dS 


(4.9) 

4i^  =  ^iCy)  =  ^Gy) =  42^ 

2/695  2/695 

(4.10) 

/v  A  (A  =  Y.V 1  (»)  -  /I  W]  Yr 

2/695 

G  (C,’2  /6 

2/e95 

=  \J2f2(z). 

Then,  it  follows  from  the  condition  uj\  <  102  that 

CIA]  =  If  /i(-A*/iK* 

z  Jds 

—  0  /  f2(-\j2f2)doj1 

1  JdS 

=  \  [  h(-*u*h)<U  +  \  [  /2(-A^/2)^2 
J  S  J  &s 

=  2  /j2(-4,/2)4 

=  7  /jW2/2|242 
j  s 

=  \  N  N  [AG)  -  f2(y)}2to2(x,y) 

xes  yes 

>  ^  N  N  lAO  -  /2(y)]Ci(a:;y) 

xes  yes 

=  \[j  V^/2|2dU 
=  44/4 

Using  Dirichlet's  principle  (Theorem  4.4)  one  sees  that  /1  =  /2  on  S. 


(ii)  In  the  proof  of  (i)  we  actually  have  proved  that  /Wl  [/1] 

f’=fi=  h  on  S 


Toj  1  [/b]  -  In  other  words,  taking 


N  N  1/4) _  /(y)]  2^2(A  y)  =  N  N  l/fr)  -  /(y)]  U  4.  y/ 

^6»S  y6»S  a;e»S  y6»S 
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or.  equivalently 

N  N  1/0*0  _  f(y)]2  •  KUy)  -  wi(x,y)]  =  0. 

xES  yES 

Therefore,  we  have 

[/0)  -  /GO] 2  •  [W2(x-.y)  -  Wi(ar,j/)]  =  0, 

for  all  x  £  S'  and  y  £  S'.  This  gives  (ii).  g 

Remark  4.3.  In  Theorem  4.1  above,  if  f  :=  fi  =  is  injective  on  S  then  we  are  able  to  get 
u>i  =  u>2  on  S  x  S.  For  example,  if  S  is  the  path  Pn  on  n  vertices  with  arbitrary  weight  ca.  then 
it  is  not  hard  to  see  that  every  nonconstant  cj-harmonic  function  /  on  Pn  is  strictly  monotonic 
and  hence  all  the  weights  are  identified.  But.  in  general,  most  graphs,  even  with  the  standard 
weight  do  not  admit  an  injective  solution  to  the  DBVP  or  NBVP.  Therefore,  it  will  be  quite 
interesting  to  figure  out  a  pair  of  graphs  and  weights  which  admits  an  injective  solution  to  the 
DBVP  or  NBVP. 

To  develop  an  idea  to  improve  Theorem  4.1  we  consider  a  graph  S  =  {1.  2.  3.  4.  5.  6}  with 
dS  =  {0.  7}  as  follows: 


* 

o 


Suppose  that  uj\  is  the  standard  weight  and  002  is  the  weight  given  by  uj\  =  002  except  only 
c^2(3.  4 )  =  k.  k  >  1.  Then  uj\  <  0J2  throughout  the  graph  S  and  uj\  =  UJ2  except  on  the  edge 
{3. 4}.  Now  define  a  function  /  :  S  — ►  R  as 

/(0)  =  a,  f(l)  =  a  -  a,  f(  2)  =  a  -  2a,  /( 3)  =  /( 4)  =  +  b)  ~  (°  +  ^ , 

/(5)  =6-2/3,  m  =6-/3,  /(7)  =  6, 

where  a.  b.  a  and  /?  are  arbitrary  real  numbers.  Then  it  is  easy  to  verify  that  /  satisfies  both 
the  equations 

A^/Or)  =  0  =  A^/O),  x  G  S. 

Here,  we  note  that  /  is  uniquely  determined  by  the  boundary  data 

/( 0)  =  a,  /( 7)  =  b, 

!U) = /<°>  -  /(i> = “■  V7) = /(7>  -  = >> 

(J(jj  I  i  Uijj  I  i 

and  each  value  f(x)  is  determined  regardless  of  the  value  c^2(3. 4)  =  k.  This  implies  that  we 
cannot  identify  the  weight  c^2(3.  4 )  =  k  even  with  all  possible  boundary  data  a.  6.  a  and  (3.  To 
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derive  a  key  idea  to  identify  LV2&  4)  =  k.  we  take  a  >  0.  b  >  0  so  that  /( 0)  >  0  and  /( 7)  >  0. 
By  a  direct  calculation  (or,  using  Corollary  2.2  )  we  see  that 

/(m)  >0,  m  —  0,  1,  2,  •  •  •  ,  7. 

Suppose  that  /  satisfies  the  relation 

(4.11)  [  fd,^  =  [  fd^  . 

Js  Js 

Then,  since 

[  fd^  =  2/(1)  +  3/(2)  +  3/(3)  +  3/(4)  +  3/(5)  +  2/(6) 

Js 

and 

[  fd^  =  2/(1)  +  3/(2)  +  (2  +  fc)/(3)  +  (2  +  fc)/(4)  +  3/(5)  +  2/(6), 

Js 

it  follows  that 

fc[/(3)  +  /(4)]  =/(3)  +  /(4), 

which  gives  k  =  1.  Therefore,  in  order  to  identify  the  weight  over  all  edges  we  need  to  impose 
an  additional  condition  such  as  (4.11). 

Now  we  return  to  the  general  situation.  We  know  that  for  a  function  :  dS  — ►  M  with 
Ids  ^  =  °  an<^  3  =  1;  2,  the  equation 


(4.12) 


'  Autjhj(x)  =  0.  x  6  S, 

<  £^(z)  =  i/>(z),  zedS., 


V  Is  h  j  ^3 

has  a  unique  pair  of  solution  (/q,  /12).  Let 

(4.13)  m.j  min 

zeds 


=  0 

hj(z),  j  =  1,2 


and 


(4.14)  mo  =  max  |my|  *  vol(S, tOj), 

where  vol(£,cOj)  =  YsxeS^'jx- 

Motivated  by  the  above  example  we  refine  Theorem  4.1  as  follows: 


Theorem  4.4.  Let  uj\  and  u>2  be  weights  with  uj\  <102  on  SxS  and  /1,  /2  :  S  — ►  R  be  functions 
satisfying  that  for  each  j  —  1,  2, 


(4.15) 


/or  a  given  function  :  dS  — » 
mo  is  the  constant  in  (4-14))- 


'A =  0,  x  6  S', 

1  SyiM  =  !«*). 

JsfjcLt  =  K 

with  fdsfj  M  0  and  a  given  constant  K  with  K  >  mo  .  (Lhere} 
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If  we  assume  that 

o 

(i)  tJi(z.y)  —  x>2 (z*,y)  on  dS  x  dS. 

(ii)  /lias  =  /2|as  > 
then  we  have 

fl=f2 

and 

ui(x,y)  =t o2{x,y) 

for  all  x  and  y  in  S. 

Proof.  We  have  already  shown  in  Theorem  4.1  that  fi=f2.  Now.  for  each  j  —  1.  2.  we 
choose  a  constant  Cj  so  that  Cj  ■  vo \{S.ujj)  =  K .  Then,  it  follows  that  Cj  >  \mj\  and.  hence. 
hj(x)  +  Cj  >  0.  x  £  S'  by  the  maximum  principle  (or.  Corollary  2.2).  Moreover,  the  function 
h{x)  hj{x)  +  Cj  satisfies  the  equation  (4.15).  By  the  uniqueness  of  the  solution  we  have 

fj(x)  =  h(x)  =  hj(x)  +  Cj  >  0.  x  G  S. 

Let  f  :=  f i  =  f2  on  S.  Then  it  follows  from  the  condition  fs  f\  dwi  —  K  —  fsf2  that 

51  f(x)dun(x)  =  N  /(^)A. 2(V 

x£S  x£S 

or.  equivalently 

51  f(x)  [V2  0)  -  (L  i  0)]  =  o. 

xES 

Since  f(x)  >  0  and  du_,]  (x)  >  dWl  (x)  for  all  x  G  .S',  we  have 

0  =  c^2  (x)  ~  doJi  ix) 

=  51  [^x-y)  -^2(x,y)]. 

eS 

Since  (x.y)  <  u>2(x.y).  we  obtain 

x>i(x,y )  =u2(x,y ) 

for  all  x  and  y  in  S.  which  is  the  required,  g 

Remark  4.5.  In  the  above  proof,  the  contribution  of  the  new  condition  that  fs  fj  d^.  —  K  >  mo 
is  used  only  but  to  guarantee  that  fj(x)  >0.  x  £  S.  Hence,  if  we  replace  this  condition  by 
/las  >  0.  j  —  1.  2  in  Theorem  4.4.  we  arrive  at  the  same  conclusion. 

As  seen  in  the  study  of  the  inverse  conductivity  problem  in  the  continuous  case  (see.  for 
instance.  [A].  [Ca],  [I],  [IP],  [SU])  it  would  be  worthwhile  to  prove  the  uniqueness  under  a 
condition  weaker  than  the  monotonicity  condition  uj\  <  uj2  imposed  the  above.  Moreover,  it 
would  also  be  interesting  to  consider  a  stability  theorem  for  the  same  conductivity  equation. 
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